The Quantum Refrigerator: The quest for absolute zero 
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Abstract 

The scaling of the optimal cooling power of a reciprocating quantum refrigerator is sought as a 
function of the cold bath temperature as T c — * 0. The working medium consists of noninteracting 
particles in a harmonic potential. Two closed-form solutions of the refrigeration cycle are analyzed, 
and compared to a numerical optimization scheme, focusing on cooling toward zero temperature. 
The optimal cycle is characterized by linear relations between the heat extracted from the cold 
bath, the energy level spacing of the working medium and the temperature. The scaling of the 

3 /2 

optimal cooling rate is found to be proportional to T c giving a dynamical interpretation to the 
third law of thermodynamics. 
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I. INTRODUCTION 



Walter Nernst stated the third law of thermodynamics as follows: "it is impossible by 
any procedure, no matter how idealized, to reduce any system to the absolute zero of tem- 
perature in a finite number of operations" This statement has been termed the 
unattainability principle 3), 0, 5), Q] . In the present study the unattainability statement is 
viewed dynamically as the vanishing of the cooling rate Q c when pumping heat from a cold 
bath whose temperature approaches absolute zero. Finding a limiting scaling law between 
the rate of cooling and temperature Qc oc T c 5 quantifies the unattainability principle. 

The second law of thermodynamics already imposes a restriction on 5 [7J. For a cyclic 
process entropy is generated only in the baths: a = —Q c /T c + Qh/Th > . If Qh 
stays bounded, \Qh\ < C, as T c approaches 0, then rearranging the inequality above gives 
C/T h > Qh/T h > Q c /T c , and so: (jr^j T c > Q c . This forces Q c — > as T c — > and, 
expanding Q c as a series near T c = 0, the dominant power 5 in Qr oc T s must satisfy 
5 > 1. Such an exponent has been realized in refrigerator models 0, [q] where the source of 
irreversibility is the heat transfer. The vanishing of Q c is also consistent with the vanishing 

n 2 k 2 T Li 

of the quantum unit of heat transport — [9]. Our goal in the present study is to set 
more stringent limits on the exponent 5 for a reciprocating four stroke cooling cycle. The 
cooling rate is replaced by the average refrigeration power 1Z C = Q c / r where r is the cycle 
period. 



II. THE QUANTUM OTTO HEAT PUMP 

We consider a refrigerator using a controllable quantum medium as its working fluid. Our 
objective is to optimize the cooling rate in the limit when the temperature T c of the cold bath 
approaches absolute zero. A necessary condition for operation is that upon contact with 
the cold bath the temperature of the working medium be lower than the bath temperature 
Tint < T c [lo|. The opposite condition exists on the hot bath. To fulfill these requirements 
the external controls modify the internal temperature by changing the energy level spacings 
of the working fluid. The control field varies between two extreme values u c and u^, where u 
is a working medium frequency induced by the external field. The working medium consists 
of an ensemble of non interacting particles in a harmonic potential. The Hamiltonian of 
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FIG. 1: A typical optimal cooling power cycle ADCB in the energy entropy Se and frequency u 
plane. 

this system, H = t^P + , is controlled by changing the curvature K = mu 2 of the 

confining potential. 

The cooling cycle consists of two heat exchange branches alternating with two adia- 
batic branches (Cf. Fig. 1). The heat exchange branches (the isochores) take place with 
uj =constant, while the adiabatic branches take place with the working medium decoupled 
from the baths. This is reminiscent of the Otto cycle in which heat is transfered to the 
working medium from the hot and cold baths under constant volume conditions. 

The heat carrying capacity of the working medium limits the amount of heat Q c which 
can be extracted from the cold bath: 

Qc = E C -E D = huj c (n c -n D ) (1) 

where Eq is the working medium internal energy at point C (Fig. [T]), Ed is the energy at 
point D and n = (N) is the expectation value of the number operator. Examining Fig. [T] 
no < rf n and no > n^, where equality is obtained under the quantum adiabatic condition 



11 



12j. This means also > n e ^, leading to Q c < fijjj c {n e c q — n e £). Maximum Q c 
is obtained for high frequency huih ^> ksTh, leading to n£ = and Ea = \tvjJh being the 
ground state energy. Then for T c — > 0: 

Q* = hu c n e c q = hu c e~^ < k B T c (2) 



where we have substituted the value of nj; 9 obtained from the partition function and the last 
inequality is obtained by optimizing with respect to uj c leading to tvjj* = k B T c . The general 
result is that as T c — > 0, Q* and u* become linear in T c . 

Only a finite cycle period r leads to a non vanishing cooling power 7Z C = Q c / r [3] . This 
cycle time r = t/, c + t c + r c h + is the sum of the times allocated to each branch Cf. Fig. 
HJ An upper bound on the cooling rate 1Z C is required to limit the exponent as T c — > 0. The 
optimal cooling rate TZ° pt depends on the time allocated to the different branches. 

The dynamics on the adiabatic segments is generated by an externally driven time de- 
pendent Hamiltonian H(u(t)). The equation of motion for an operator O of the working 
medium is: 

Typically [H(t), H(i')] ^ which leads to friction like phenomena [14|,ll5j]- Too fast adiabatic 
segments will generate parasitic internal energy which will have to be dissipated to the heat 
baths, thus limiting the performance. The external power of the compression/expansion 
segments is the rate of change of the internal energy of the working medium [if]]. Therefore 
inserting H for O in Eq. ([3]) leads to the power ^| = V = (^jr)- The dynamics on the 
adiabatic segments is unitary, therefore the von Neumann entropy S vn = —k,Btr{p\n.p) 
is constant. In contrast the energy entropy Se changes, where Se = — fcs Pj In .P/ and 
Pj = tr{\j)(j\p} is the probability of occupying the energy level j. Constant Se is obtained 
only under quasistatic conditions. 

Using the Heisenberg picture, the dynamics on the heat exchange branches, termed iso- 
chores, are generated by £*(0) = |[H, 6] + C* D (0) [l7] with the dissipative Lindblad term 
C* D leading the system toward thermal equilibrium of an harmonic oscillator defined by 
jr = exp(— -j^f) l^L- Fo r the dissipative dynamics, the heat flow from the cold/hot bath is 



Q = <£d(h)> HJ,H 

At thermal equilibrium the energy expectation value is sufficient to fully characterize the 
state of a system. For the working medium not in equilibrium, there is a family of generalized 
Gibbs states [jj]] that completely characterize the system during the cycle. These Gibbs 

1 * 2 K(t) ~ 2 

states are defined by three operators: the time dependent Hamiltonian H = ^-P H — ^-Q , 
the Lagrangian L = — -^p-Q and the correlation C = cj(t)|(QP + PQ). As a result 

p = p(H,L, C). Starting from any initial state, the system will reach a generalized Gibbs 
state and will remain in such states. The invariance of the set H, L, C under the equation 



4 



of motion, is due to this set forming a closed Lie algebra, which leads to closed equations of 
motion on the adiabats as well as on the isochores [la . Il9| ] . 

The dynamics of the operators on the adiabats is obtained from Eq. ([3]): 
/tt\ / „ _.. n \ /tt\ 

d 



dt 
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where ji = is the dimensionless adiabatic parameter. The power becomes: P = //u;((H) 



L)) [15|, [l9|. The solution of Eq. (jlj) depends on the functional form of u;(£). When 



/i < 1, the number n(£) will remain constant on the adiabats; these are the quasistatic 
conditions. For most other functions ou(t), the time evolution will involve some quantum 

n 

friction [15J and n/ > rii due to the resultant parasitic increase in the internal energy 
AE = hujf(nf —Hi). The dissipation of this energy in particular into the cold bath counters 
the cooling: Q c < hw c (nl q — n c ), therefore when n c > rf c q the refrigerator can no longer cool. 
On the isochores the energy displays an exponential approach to equilibrium: 

— = -r(H-(H) eg i) (5) 

where r = k\ — fcf is the heat conductance. (H) eg is the equilibrium expectation of the 
energy. The heat transfer becomes: Q = — r((H) — (H) eg ). 

The operators L and C display an oscillatory decay to an expectation value of zero at 
equilibrium: 

S(')«=(L r -r)tS 

The equation of motion (j4j), (JHJ) and (EI) can be solved in closed form for certain special 
choices of u(t) (Cf. section [mi below) and numerically for any given functions u(t) and time 
allocation to the branches. After a few cycles, the refrigerator settles down to a periodic limit 
cycle which allows to calculate the cooling power TZ C = Q c /t from the expectations of 
H, L, C in the limit cycle. 



(6) 



III. OPTIMIZATION OF THE COOLING RATE 



For sufficiently low T c , the rate limiting branch of our cycle is cooling the working medium 
to a temperature below T c (A— >D along the expansion adiabat). As T c — > 0, the total cycle 



U (t) - w ® 



time t is of the order of the time of this cooling adiabat, Th c , which tends to infinity. 

Quantum friction is completely eliminated if the adiabat proceeds quasistatically with 
// < 1. This leads to a scaling law 1Z C oc T s with 5 > 3. It turns out however that it is not 
the only frictionless way to reach the final state at energy Ed = (w c /wft)i?A- We describe 
two other possibilities which require less time and result in improved scaling, 5 = 2 and 
5 = 3/2 respectively. 

The first frictionless solution to Eq. (jlj) is obtained for fi = const, by changing the time 
variable to 9 = j^uo{t')dt'. Then factoring out the term fxl and diagonalizing the time 
independent part with the eigenvalues Ao = and \± = ±Q where Q = a//x 2 — 4 leads to 
the adiabatic propagator U a of H, L, C: 

/ /i 2 c - 4 fxQs 2/i(c - 1) ^ 

fxQs tt 2 c 2fis • (7) 
\ -2/i(c - 1) -2Qs fi 2 - 4c , 

where c = cosh(f26 l ), s = smh(fl9) and 9{t) = — log(^^)//i. The cycle propagator becomes 
the product of the segment propagators U cyc = U. c U c hUbUhci where Uh/ C is obtained from Eq. 
([5]) and Eq. ([6]) on the isochores. 

The energy change on the expansion adiabat is the key for the optimal solution: A — > D 

E D = l -hw c cosh(fi0 c ) - 4) ,B C =-- log (C) (8) 

where C = ^ is the compression ratio and equilibration is assumed at the end of the hot 
isochore Ea = \hjJh for Uh — > oo. For very fast expansion \i — > oo, Ed = jHuj c (1/C + C). 
As T c — > 0, Ed = \fajJh which becomes larger than E* q therefore the cooling stops due to 
friction. For the limit of infinite time /i — > leading to the frictionless result characterized 
by constant n and Se- Then Ed —> \^> c which is the ground state of the oscillator. At this 
limit since r — >• oo, 1Z C = 0. The surprising point is that we can find an additional frictionless 
point where n c = nt, when cosh(Q9 c ) = 1. Then /i < 2 and Q becomes imaginary leading 
to the critical points: 

2 log (C) 

fi = ===== (9) 



'47T 2 + log (C) 2 

and t£ c = (1 — C)/ (n*Uh)- Asymptotically as T c — > and u; c — > 0, the critical terms approach 
/j,* — > —2 and with it the time allocation r^ c = I^J 1 . This frictionless solution with a 



minimum time allocation r^ c scales as the inverse frequency uo^ 1 which is better than the 
quasistatic limit where r hc oc oj~ 2 . As we will see, it leads to 5 = 2. 

Inspired by these findings, we sought the minimum time frictionless solution. The result- 
ing optimal control problem [^(J is solvable leading to a second closed form solution. The 
optimal trajectory is of the bang-bang form with three jumps 



oj h , for t = 
lu c , for < t < Ti 

(10) 

u h , for n < t < r hc 



for t = r h 



C 1 



where n + r 2 = r hc and the times n = ^ arccos ( ( "f^f) 2 ) and r 2 = ^ arccos 

are chosen such that the number operator is preserved n/ = n^. The minimum time allo- 

_i 

cation for uj c — > which is appropriate for T c — > becomes rj* = -t^u; c 2 , which is better 
than the solution in Eq. (??). As we show below, it leads to 5 = 3/2. 

Both frictionless solutions lead to an upper bound on the optimal cooling rate of the 
form: 

K c < Au v n e c q (11) 

where A is a constant and the exponent v is either v — 2 for the /i = const solution or v — | 
for the three-jump solution. Optimizing TZ C with respect to uj c leads to a linear relation 
between u c and T c , ftw c = KkgTc. The constant k = 2 + 7 ? (— 2e~ 2 ) ~ 1.6 for v = 2 and 
k = 3/2 + "P(— 3/2e~ 3 / 2 ) « 0.87 for z/ = |, where "P is the product-log function. 

Once the time allocation on the adiabats is set the time allocation on the isochores is 



optimized using the method of Ref . [15] : 



K = J^y^cK" - n h q ) (12) 

where z = ThTh = T c r c and z is determined by the equation 2z + r(r/ lc + r^) = 2sinh(z). 
For the limit T c — > 0, Tthc is large therefore z is large leading to: 

K - ^g^ rto c ^ - n?) (13) 

(1 + 1 [J hc + r c/l J) z 

At high compression ratio Uh 3> ^ c and if in addition w c < T we obtain: 

K « (14) 



for the [i = const frictionless solution, and 

K* c « Uwl^n^ , (15) 

for the three-jump frictionless solution. Due to the linear relation between u c and T c , Eq. 
( I14p and (fT5|) the exponent 5 where 5 = 3 for the quasistatic scheduling, 5 = 2 for the 
constant \i frictionless scheduling and S = | for the three-jump frictionless scheduling. 

To check the optimization assumptions a numerical procedure was applied to maximize 
the cooling rate by adjusting the times on the four branches for a given a choice of scheduling 
function and the external constraints on the cycle. These constraints are the coupling T, 
the temperatures T c and T h , and the frequencies u c and ojh- The cooling rate optimizations 
employed random time allocations to the different cycle segments augmented by a guided- 
search algorithm. The choice of scheduling function u(t) determines the exponent of the 
scaling in TZ C oc T s . The optimal cooling rate for linear and exponential scheduling functions 
are shown in Fig. As a final numerical corroboration, we tried a multistep genetic 
algorithm allowing piecewise variation of u(t). The algorithm converged to a cooling rate 
very close to the optimal three jump solution. 

Two main observations have led to the optimal exponents as T c —>■ 0, the first is that 
the time allocation on the expansion adiabat sets the scaling and the second is that the 
frictionless cycles have superior performance. Fig. [2] also shows the results of numerical 
optimizations for the two frictionless schedules. At low temperatures the time allocated to 
the adiabats dominates and scales as r^ c oc 1/T C for the \i = const, schedule and r^ c oc 1/T<V 2 
for the three-jump schedule. Since Q c for all cases is linear with T c , the asymptotic cooling 
rate approaches 1Z C oc T c 2 and TZ C oc T c respectively. 



IV. DISCUSSION AND CONCLUSION 



The optimal quantum refrigerator in the quest to reach the absolute zero temperature 
shows a linear scaling of Q* with u c and T c . This scaling is the minimum to eliminate the 
divergence of the entropy generated on the cold bath. If the energy level spacing hu c cannot 



follow T c the refrigerator will be limited by a minimum temperature [2l|. If the level spacing 
follows T c , the scaling of the cycle time is dominated by the scheduling function u(t) on the 
adiabats. The best results were obtained for the three-jump frictionless solutions which give 
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FIG. 2: The cooling rate as a function of temperature for different scheduling functions u>(t). 
Straight lines are linear continuations of the last points. The exponent of 1Z C oc T% is indicated. 
The lowest exponent (red squares) shows the three-jump frictionless optimization. The next lowest 
(blue circles) show the result of the optimization with [i = const. The black diamonds correspond 
to u(t) oc exp(at) and the highest exponent corresponds to u(t) oc t (magenta triangles). 

t oc co> c -1 / 2 . The three-jump scheduling is the minimum time frictionless solution [20]. We 
conjecture that any cooling cycle is limited by the adiabatic expansion [5]. Our conjecture 
implies that the unattainability principle is a consequence of dynamical considerations and 
is limited by the exponent lZ° c pt oc T c 3//2 . 
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